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Abstract. Dark Matter (DM) theories and mass-tracing-light theories like MOND are by construction nearly 
degenerate on galactic scales, but not when it comes to the predicted shapes of Roche Lobes of a two-body 
system (e.g., a globular cluster orbiting a host galaxy). We show that the flattening of the Roche lobe is sensitive 
to the function fj,(g) in modification of the law of gravity. We generalise the analytical results obtained in the 
deep-MOND limit by Zhao (2005, astro-ph/0511713 ), and consider a binary in the framework of a MOND-like 
gravity modification function /i(g) or a general non-Keplerian gravity g oc . We give analytical expressions 
for the inner Lagrange point and Robe lobe axis ratios. Fhe Roche lobe volume is proven to scale linearly with 
the true mass ratio, which applies to any g.(g), hence mass-tracing light models would overpredict the Roche 
lobe of a DM-poor globular cluster in a DM-rich host galaxy, and underpredict the size of a DM-richer dwarf 
satellite. The lobes are squashed with the flattening ~ 0.4 in the strong gravity and ~ 0.6 in the weak gravity; 
a precise measurement of the flattening could be used to verify the anisotropic dilation effect which is generic 
to MOND-like gravity. We generalise these results for extended mass distribution, and compare predicted Roche 
radii with limiting radii of observed globular clusters and dwarf galaxy satellites. 
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^ ' 1. Introduction 

h : 

5t i The law of gravity is still uncertain experimentally on very small scales and very large scales. An example is the 
running controversy (Sellwood & Kowsowsky 2001, 2002 and references therein) over dark matter or MOND (Modified 
Newtonian Dynamics) as the right explanation of why the Newtonian gravity g^ from baryonic stars and gas in galaxies 

fall short of explaining the observed acceleration g by a factor fj,(g) — ^ a ™ +gn ) " ~ (0.01 — 1), where ao ~ 10~ 8 cmsec -2 

is a characteristic accerleration, and n = 2 (see examples in Sanders & McGaugh 2002) or n — 1 (Zhao & Famaey 2006). 
This is normally the justification to invoke dark matter particles on galaxy scales so that M/L ~ 1/fi ~ (1 — 100); the 
mass profile of the dark particles needs to be well-coupled to the baryon distribution (McGaugh 2005). Alternatively 
often a good fit can be made in a MOND gravity (Milgrom 1983) or MOND-like gravities (Bekenstein 2004, Sanders 
2006). The important question for astronomers is whether these two views of fundamental difference are by and large 
equally good descriptions of features in galactic systems (e.g., the orbits of halo streams in the Milky Way, Read 
& Moore 2005), hence are degenerate, or there are structures in galaxies which simply cannot be explained by any 
modifications of the law of gravity (e.g., problems on sub-galactic scale, Zhao 2005). 

In MOND the gravity g = — V$ and gravitational potential <& follow a Poisson-like equation 
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£M,*(Ri), Sfi = sK\s\), (1) 



for an isolated system of N baryonic point masses Mi, M/v at distances i?,;. In general, the vector g^ is not curl-free, 
and equals the curl-free Newtonian gravity gw plus a divergence-free curl-field. Here < fi(g) < 1 is the modification 
factor, and is a function of the amplitude of the gravity g = |g| (or equivalently an implicit function of g^ = | g M | ) . 
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if fi = 1, then g = g^ = gjy and the normal Poisson's equation is recovered. This modification satisfies the usual 
conservations of total energy and angular momentum for any geometry of baryonic distribution (see Bekenstein & 
Milgrom 1984). 

In the MOND theory and in spherical symmetry, the gravity g(R) at distance R from a baryonic point mass M is 
boosted from the Newtonian value gN(R) = GM/R 2 by a factor > 1 such that far away from the mass point we 
have 



9 ju 9n/o,q ( g \ ( R 



g/a \aoJ \ Ro 



(2) 



This reproduces the flat rotation curves in bright disk galaxies at radii beyond i?o = J since V 2 ~ gR ~ aoRo 
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est. The predictive power of this 20- year-old classical theory with virtually no free parameters (Bekenstein & Milgrom 
1984) is recently highlighted by the astonishingly good fits to contemporary kinematic data of a wide variety of high 
and low surface brightness spiral and elliptical galaxies; even the fine details of the ups and downs of velocity curves 
are elegantly reproduced without fine tuning of the baryonic model (Sanders & McGaugh 2002, Milgrom & Sanders 
2003). For a long time a problem has been that the gravitational lensing in this non-relativistic theory is illly-posed, 
and it misses a factor of two in the bending angle (e.g., Qin, Wu & Zou 1995). However, this empirical MOND has 
now a respectable relativistic field theory formulation (called TeVeS by Bekenstein 2004), which passes standard tests 
to check General Relativity, and allows for rigourous modeling of Hubble expansion and gravitational lensing. This 
has generated wide interests, and many are examining the consequences of modifications to gravity (Skordis, Mota, 
Ferreira et al. 2005, Hao & Akhoury 2005, Chiu, Ko & Tian 2005, Zhao, Bacon, Taylor et al. 2006, Pointecouteau & 
Silk 2005, Ciotti & Binney 2004, Ciotti, Londrillo, Nipoti, 2006, Baumgardt, Grebel, Kroupa 2005, Read & Moore 
2005, Famaey & Binney 2005, Zhao & Famaey 2006). 

Zhao (2005) has shown that the shape of the Roche lobes in the deep-MOND limit is more squashed than in the 
Newtonian case. Here we study the shape of the Roche lobe of a two-body system in a MOND-like gravity. We show 
that Roche lobes shape varies with the assumed law of gravity, and could be used to differentiate among laws of 
gravity. The results are first shown for point masses and are generalised to extended mass distribution. 

2. MOND-like gravity 

A very broad class of modified gravity models could be parametrized by a modification function given as follows, 
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(3) 



which makes /i(g) an implicit function of the gravitational strength g = |g|. The Newtonian gravity corresponds to 
models with k = n = 0. Exotic gravity can be achieved by letting, e.g., (k, n) = (1, 1) 1 or (k, n) = (— f , l) 2 - Example 
rotation curves are shown in Fig. ^ for a point mass (upper left panel) and an extended mass (upper right panel) . 
Conventional MOND gravity corresponds to models with k = 1/2. Such a model in spherical symmetry goes from 
9 ~ 9n ~ 9n = GMR~ 2 in strong gravity to g ~ ^/gWao = \ZGMoqR~ 1 in weak gravity around a point mass M. The 
sharpness of the transition is controlled by the parameter n. The location of the transition zone for a point mass M 

j~GM 

can be inverted as 

-2 In 

i / i 

fi(x) = x 



is given by Rq = J To see the link with conventional MOND fj, function, we note that for k — 1/2 equation J3J 



x = — , for k = -, (4) 
a 2 



which approaches i or 1 for small or big x. Our function in the case (k, n) = (1/2,3) approximates the standard 
MOND function p = x(l + x 2 )^ 1 / 2 very well (cf. bottom panel of Fig. [T}. And the case with (k,n) = (1/2,3/2) 
approximates the simple MOND function fi = x(l + x)^ 1 recommended by Zhao & Famaey (2006). 

For more compact notations later on we shall introduce a few auxilary functions, related to fi in a spherical 
potential. We define 

A ^) = 77tA A 1 = A + 1 = — (5) 
a In g a In g 

1 so that the gravity g — > est at large distances as in, e.g., (Mannheim 1997). 

2 so that g cx at small distances as in, e.g., the theory of Qin et al. 2005 with three extra nanometer-thick spatial 
dimensions. 
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Fig. 1. Top panels: 
(cf. eq.EJ) and V cir (R) 



x=g/aO 

shows the circular velocity curves V c i r (R) = \fg~R in various modified gravity (red circle) 
= \Z9nR in Newtonian gravity (blue cross); both in units of (GMqcio) 1 ^ 4 . These are shown 
in left panels for point mass with (k,n) = (1/2,1) (smallest symbols), {k,n) — (1/2,8) (bigger symbols) and in 
right panels for an extended Hernquist profile mass with a scale b — 1/2 in a stronger-than-MOND gravity with 
(k,n) = (9/10, 1). Middle panels: Same parameters as in upper panels, except showing £ (red circles) and A (black 
diamonds) as functions of rescaled orbital distance R/Rq. Bottom panel: compares the function n(g) in the popular 
MOND models (/J,(x) — x/y/l + x 2 in boxes, and x/(l + x) in circles) with our function in the case (k,n) = (1/2,3) 
(thick line) and (1/2,3/2) (thin line) and (1, 1) (dashed line). 



We also define 

'dlnR(g) 



CO?) 



G?ln g 



Ci ^ C + 1 
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where fl(R) = ^/g is the angular frequency for a circular orbit of radius R. In writing down R(g) we treat the spherical 
radius R as a function of the gravitational field strength g through the relation MR~ 2 oc g^ ~ g^ = [ig. Note in the 
case of the field around a point mass, and assume non-negative n and k, A(g) is related to the parameter £(g) by 



ding 
d\nR-- 



1 

AT 



1 near a point mass 

1 — k far from a point mass. 



(7) 
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Conventional MOND models have the property that (A, £) — > (0,2) when the gravity is strong and (A, £) — > (1, 1) 
when the gravity is much weaker than clq. The result in this paper applies to any function A(g) or £(<?)• In general, 
we shall call the above gravity g the MOND-like Gravity or Modified Gravity. 

Fig- HI shows how £ (red circles) and A (black diamonds) change as the MONDian gravity decreases away from a 
point mass (middle left panel) and from an extended mass (middle right panel). These MOND models have a narrow 
range of 1 < C < 2 or 1 > A > 0. 



3. Roche Lobes of a two-body system in a MOND-like gravity 

A binary system is specified by two baryonic masses M and m separated by a distance D a along, say, the z-axis. The 
two bodies could be a galaxy and a satellite. The MOND-like gravity distribution is determined by equation ^ where 
N = 2 for a binary. Define the inner Lagrange point of the binary at a distance r\ and D D — n from the masses m 
and M. There are only 3 independent dimensionless quantities in this system, r\/{D — r\), m/M and a typical value 
for the parameter £. From dimensional arguments we expect that the Roche radius might be expressable in general as 



t i I m x 



D \A(C)M 



(8) 



where A and B are dimensionless functions with dimensionless arguments; we now go on to show, using perturbation 
analysis, that this is indeed the case, and we set the exact expressions for A and B. 

We set up the coordinate system such that the origin is on the smaller mass m, at a distance D from the mass M 
on the z-axis. Let the low- mass satellite with m/M <C 1 rotate around the mass M (fixed) with an angular velocity 
fly, then particles in the corotating frame conserve the Jacobi energy with an effective (triaxial) potential 

<t> e (x,y,z) = $- x2 + {z ~ Do)2 n 2 , $ = * + *i, (9) 

where the gravitational potential <£> of the galaxy-satellite system is expressed as two terms, where we treat the effect 
of the smaller mass m as a perturbation $i to the unperturbed potential $0 of the bigger mass M. 

First consider the unperturbed case where we set the satellite m — 0. At a distance R far way from the isolated 
baryonic mass M, the gravitational potential can in general be approximated as spherical with 



$0 = / g(R)dR, R = yJ(z-D )2 + y 2 + x* (10) 



where the gravity g(R) at radius R from the spherical galaxy centre is related to the Newtonian gravity — GMR/il 3 

by 



S^ = SN = - — [x,y,z~D } 1 g = -V*o(-R)- (11) 



GM 

Taylor-expand the three components of the gravity in the vicinity of the origin (x, y, z) = (0, 0, 0), we have 



v$ ^,y,A^ 2 + [o,o,z]^^, (12) 

The gravity field is such that if a test particle is put on circular orbit with an orbital radius R = D around M, 
then the rotation angular frequency must be Q in order to balance gravity with centrifugal force. Integrate the above 
potential gradient, we find 

$ (i?) w $ (^ o ) + [x 2 + y 2 - z 2 (] ?L + zD o n 2 (13) 

in the vicinity of the origin (0, 0, 0), where according to the definitions of ( and A we have £ = — d |j°o 3 ^'°' > ■ 

Now we consider the perturbation due to the satellite. Milgrom (1986) showed that in a medium under a uniform 
external field g^z with a dielectric index [i(go), the perturbation due to the added point-mass m is almost Newtonian 
apart from a mild anisotropy with a potential (cf. his eq. 11) 

G 



$i(x,y,z) = — , f= v/z 2 + Ai(y 2 + x 2 ), (14) 

where m! — mj '/it (go) i s the modified inertia of the point mass, and f is the effective distance from the centre of the 
satellite and Ai = A + 1 = 1+ dl °s^(g°) ^ l is the MOND dilation factor. Here we show that we can recover the same 

1 d log SO ' 

expression for a satellite in a galaxy. Sufficiently far away from the satellite, the gravity is dominated by the spatially 



HongSheng Zhao and LanLan Tian: Roche Lobe Shapes for testing MOND-like Modified Gravities 




Fig. 2. shows one octant of the 3-dimensional Roche Lobe of a satellite of mass m rotating in the xz plane around 

/ \l/3 

a mass M at a distance D Q on the z-axis. The sizes are scaled by a length s = ( M ^ ) J D Q . The curves are for 
C = - f, '"n = 2 (Newtonian, thick black) and C = 1 (Deep-MOND, red). A slimer Roche lobe for £ = 1/5 (thin green) 



and a fatter Roche lobe for £ = 5 (thin blue) are also made for comparison. 



slow- varying gravity field of the unperturbed galaxy, which we denote with the shorthand gR = |g(R)|. We also use 
the shorthand /ir = ^(|g(R)|) to denote the slow-varying MOND "dielectric index" due to the galaxy. The perturbed 
dielectric index is /xr plus a perturbation due to the self-gravity of the satellite, which is given by 

M(|g|)-MR«^[|g|-PR]«^R-V$i. (15) 
ag ag 

The perturbation part of the modified Poisson's equation (eq. ^) becomes 

V ™ V (^ R+ ^ VR$1 ) (gR ' V<&l) ., XM fl< n 

V • _ - V — « mS(r). (16) 

Cancel out the Oth order terms due to the mass M only, drop terms of order |V$i| 2 and small terms of order -^V$i, 
we end up with a linear equation for $i, which in the vicinity of the satellite is given by 

(Aid 2 + d 2 x + d 2 y ) $i = 4wGm'd{r), (17) 
where 



Gm' 1 g{D ) _ n 2 D, 



Gm ^ R g^Do) GMDq 2 



(18) 



Apply the transformation z' = z\J Ai, x' = x, y' = y, it can be shown that to the first order in m the perturbation in 
potential is given by equation (|14fl . where go f=a $r ~ g(Do). As verified by Milgrom (1986) using a numerical solver, 
our equation (14) is "a good approximation" as long as the Newtonian force of the satellite is small compared to 
the Newtonian part of the external field gofJ-(go): "thus the equipotential surfaces are ellipsoids of revolution with the 
major axis in the direction of the external field." 

The combined potential <f>(x,y,z) — <&o{R) + ^i{x,y,z) clearly has axial symmetry around the z-axis, and has 
two centres separated by distance D a along the z-axis. The above formulation approximates the potential of, e.g., the 
Milky Way galaxy of baryonic mass M plus a satellite of baryonic mass m. 3 



3 To check the validity of our approximations, we also compute the density p from the modified Poisson's 
equation p = — V • g M , and the mass m(r) enclosed in a radius r from G ™>f = ~ J Tw ' Sc* ^ e nn d that the 
density p ~ lO" 10 -^ ~ 0, and - 1 ~ iS« ~ unless very close to the smaller pointmass m ~ (1CT 6 - 10~ 3 )M 
or the bigger pointmass M. The approximations are less good if the bigger mass M is extended. 
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The effective potential <fr e {x, y, z), if Taylor expanded near the secondary (0, 0, 0), has an triaxial shape with 
y 2 - (iz 2 D m/M 



$ e w est + Cl z Di 



mi 



^z 2 + A 1 (y 2 +x 2 ) 



(19) 



where the terms to the first order of z and second order to x are canceled due to force balance. The inner or outer 
Lagrangian points is then calculated from the saddle point of the effective potential where 







d$ e (0,0,z) 



dz 



= ci 2 d 



riC 



D Mr 2 



This definition of the Lagrange radius n reduces to 
M 



(20) 



(21) 



So inside the Lagrange radius T\ the average density of the satellite equals k-times the average density inside the 
orbital radius D a of the satellite. The equation can also be written as 



n 

Do 



CiM 



Ci = i + C- 



(22) 



Note that the masses here m and M are true baryonic masses of the binary, not the modified inertia masses. Eq. I|21|l 
can be further massaged into 



IGt 



(23) 



where CI' takes the meaning of the angular frequency of circulation around the secondary at its Robe radius. Hence 
the tidal radius of the secondary scales linearly with the terminal velocity of the secondary, and the internal period is 
within a factor of of the period of the secondary's orbit. Finally note that near the tidal radius, the gravity due 
to the secondary and primary has the ratio 

Gm' / m\i , , 

ims-^*) (24) 

which confirms that near the tidal radius n it is valid to approximate the self-gravity of a low-mass secondary 
(m/M <1) as a perturbation to the slow- varying gravity g(D a ) of the primary. 

The shape of the Roche Lobe is specified by contour of equal effective potential with the contour passing the 
Lagrange points (0,0, ±ri). The effective potential can be redefined as a dimensionless function ty e (x,y,z) after 
massaging zero point and a prefactor. Substitute (m/M) in favor of r\/D using eq. I|22|) . we found that the points 
(x, y, z) along the Roche lobe satisfy 



*e = 



3Cxr| + y 2 - Ci* 2 
2r 2 



Cm 



V^ 2 



Air 



= 0. 



(25) 



Clearly the z-axis intersects the Roche Lobe at (0,0, ±ri), where 

n / m 
D~ = \(^M 

which defines the first axis of the Roche Lobe. Draw a line parallel to the x-axis direction, it intersects the Roche lobe 
at (±r 2 ,0,0). Solving equation (|2*5|l we find 

3r 2 1 

2n 



(26) 



(27) 



/AT' 

which defines the second axis of the Roche Lobe. Likewise the third axis of the Roche Lobe is defined by the intersection 
points (0, ±r3, 0) with a line parallel to the rotational y-axis. We find 



-i = ( U a -l* tl + ll-tt-l! 

3r 2 L 



\A + AxCi 



by simply solving equation (|25|l . which reduces to an essentially cubic equation 

2 



— - 1 = 

r 3 



27AiCi 



(28) 



(29) 
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-dLog g/dLog R 

Fig. 3. shows the Roche lobe semi-axis sizes ri,r2,r3 of a satellite m at orbital distance R = D as functions of 

1/3 

the gravity power-law index £ = — ^ogD • Here the sizes are rescaled by s = ^ m{d ) ) Also shown are the 

arithematic mean (ri +T2 + t"3)/3 (thicker dashed curve), and the geometric mean size l^i^rs! 1 / 3 (thin dashed curve). 

4. Generalisation and Discussion 

4.1. Extended spherical primary mass 

Although we refer to our primary mass as a mass point, our result applies exactly for a spatially extended spherical 
mass distribution of the primary M(R), e.g., a galaxy, if we simply make the substitution 



M -> M(R). 

As a result the ( — A relation changes to 

d\nM . 
2-- r ^=CA 1 . 



d\nR 

As a specific example, a host galaxy can be approximated by a spherical Hernquist mass profile with 

2 



M(R) = M 



R 



R + VhRc 



R<\ 



GM 
a 



(30) 
(31) 

(32) 



with the scalelength parameter yh — 1/2, where we adopted notations of Zhao, Bacon, Taylor et al. (2005). The 
Newtonian gravity of the model is 



GM(R) 2 
9n = — ^ — = a ° \V h + V) 



R 

Ho 



(33) 



Note gN — »o/(l + b) 2 — 4ao/9 < ao at R — Rq and = o-o/b 2 = 4ao is finite at R = 0. We give the Hernquist 
model a stronger-than-MOND gravity with (k,n) — (9/10, 1) so that the model has a rising terminal velocity curve. 
As shown by Fig. lb this model has £ = near the primary, and £ = 1/5 far away. 



4.2. Prolate or oblate? 

Interestingly, the second and the third (i.e., rotation) axes of the Roche lobe are always very similar in length. This is 
because equation l|29|) contains a factor 27^^ = WTKTSs ^ 1 f° r an y — 1 < A < 00. In fact that 



27(A+3) 



[1, 1.02, 1.03, 1.045, 1.055, 1.065] w 1, 



(34) 
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Fig. 4. similar to the previous one, but shows the aspect ratio of the Roche lobe semi-axis sizes rijr\ , r^/ri as 
functions fC = -|g* 





Fig. 5. showing the Roche lobe axis ratios ri/ri, r^/ri as functions of R/Rq with R = D being the orbital distance 
of the secondary. Left panel: Models are as in Fig. la and Fig.lc, with (k,n) = (1/2,1) (solid) and (k,n) = (1/2,8) 
(dashed). Right panel: as in Fig. lb and Fig. Id, with (k,n) = (9/10, 1) (solid) and (k,n) — (1/4, 1) (dashed) around 
a Hernquist profile mass. 



for £ = [0, 0.5, 1, 2, 4, oo], where (=1,2 are the deep-MOND, and deep-Newtonian cases. 

The first axis is often substantially longer than the other two, but not always. We find for point masses 

— = [oo, 3, 2.12, 1.5, 0.95,0], for C = [0,0.5,1,2,5,0c]. 



(35) 



The generally triaxial Roche lobe for point masses is nearly prolate for < ( < 5, and it becomes nearly spherical 
for £ = 4 — 6, and becomes nearly oblate only in exotic situations where the power-law index £ = \d\ng/d\riR\ is 
"harder" than 5. For extended mass, the ratio can be calculated from 



n 

r-2 



3 / rilogffjy 
2V dlogg ' 



(36) 



together with a relation between g and g^ (cf. eq. EJ) and g^ — G ^i R ^ ■ 

Overall the Roche lobe axis ratio changes by a significant amount for different assumptions of the laws of gravity 
and the mass profile (cf. Fig. [SJ, e.g., the MONDian Roche lobe (k = 1/2) at large distance is significantly more 
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Fig. 6. showing the Roche lobe semi-axis sizes n, ^Jr^r-i, as functions of R/Rq. Here the sizes are rescaled by s = 
( \ 1/3 

( m\h) ) R W1 th R = D being the orbital distance of the secondary. Left panel: models as in Fig. la and Fig.lc, 
with (k,n) = (1/2,1) (solid) and (k,n) — (1/2,8) (dashed). Right panel: models as in Fig. lb and Fig. Id, with 
(k, n) = (9/10,1) (solid) and (k,n) = (1/4,1) (dashed) around a Hernquist profile mass. 



squashed than the Newtonian case, and a modified gravity with k = 9/10 at large distance is even more squashed. 
We speculate that it might be feasible to differentiate these laws with a precise measurement of a Roche lobe. One 
advantage of using axis ratio is that it does not require precise knowledge of the satellite baryonic mass and its 
actual distance and absolute size; none of these factors enter equations (|27|l and (|28l) . In particular we note that in 
Newtonian gravity Ai = 1, hence = | holds for any mass distribution at any distance; it would be interesting for 
testing non-Newtonian gravity to quantify the deviation from this ratio in observed Roche lobe. 



4.3. Size-mass and size-distance relations 

In general, the Roche radii of two-body systems (whether stellar binaries or satellite-host galactic binaries) provide a 
measure of the relative weights of the binary in baryons in any MOND-like gravity. It is interesting that the Roche 
Lobe size in any MOND-like gravity has the same simple scaling relation 

m , M(D ) _ dlogVl 

^ /- ^T ~ Cl ^glT (37) 

The above relation holds for an extended primary as well. As in the Newtonian case the relevant mass is the unmodified 
true baryonic mass of the binary, the difference is only in the prefactor, which depends on the M and D . For point 
masses the prefactor varies from d = 3 in the region with Keplerian rotation V^; r oc Rr 1 !" 1 (Binney & Tremaine 1987) 
to (i = 2 for the flat rotation curve region. A similar relation holds for r2, 

-1/3 

OlR = 



D \MJ 3A 2£3 



where 



t d\ogV c l dlogg N 
a log it d log g 

Here ©ir takes the meaning of a rescaled Roche lobe intermediate semi-axis. 

The size of the lobe varies with the distance from the primary. The Roche Lobe sizes in terms of n and y/rWE ~ f 2 
can be read off from Fig. Elfor different assumptions of the /x-function and mass profile; in practise we can treat the 
Roche lobe as prolate. 



5. An application: Limiting radii of outer satellites of MW and M31 

As an illustration of the usuage of the Roche lobe, we compare the predicted Roche sizes with the observed limiting 
radii of outer satellites of galaxies. Slightly different from Zhao (2005), the sample consists of all globular clusters and 
dwarf galaxies beyond 20 kpc (outside the disks) of both the Milky Way and Andromeda. Only Andromeda satellites 
with reliable 3-dimensional distances to the Andromeda centre are used (McConnachie et al. 2004); tidal radius and 
distance of the newly discovered AndIX are from Harbeck et al. (2005). The limiting angular radii #ii m! mi n are taken 
primarily from Harris (1996) and Mateo (1998). Most systems are consistent with no flattening in projection. For a 
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few significant flattened systems we take the limiting radii on the minor axis. This radius should approximate the 
intermediate semi-axis r-i of a nearly prolate Roche lobe if Roche lobes are filled; a prolate object appears the same 
width on the minor axis at any inclination angle. Hence we have 

7T = ^ 1/3 ^^"' Fx<l, (40) 

L>o JJq JJ q 

where the factor F\ < 1 takes into the facts that (i) the Roche volume is smaller at the pericenter of an eccentric 
orbit of a satellite, and that (ii) the volume filling factor of the pericentric Roche lobe can still be less than unity if 
star counts drop below detection limit before reaching the edge of the lobe (e.g., Grebel et al. 2000). Several distances 
are involved here, D is the present orbital distance, D is the present distance of the Sun from the satellite, and Dq 
is the Sun's distance to the centre of the host galaxy (here the Milky Way or the Andromeda). The rhs is insensitive 
to errors on distance D; for the outer Milky Way satellites we have -jj- ~ 1, and for Andromeda satellites we have 
■p^ ~ 1. Assuming mass-tracing-light, we expect that the satellite-to- host mass ratio is the ratio of observable stars 
within a factor of two, allowing for some differences in age and metallicity. So we can define 

Tfl 

F 2 = — 0.5 - 2. (41) 

Multiplying a factor (m/M)^ 1 / 3 to both sides of the above equation (|40(1 . we have 



''2 

D 



^)~ V W 1/3 e^<e roin , (42) 



where 

JJq JJ 

is an observable, taking the meaning of the rescaled limiting minor axis size of a satellite. If MOND-like mass-tracing- 
light models are correct, substitute in equations (|26|l and l|27|) we expect that 

2 2 1 / 6 

e min > m — ~ — = 0.374, (44) 

3 Aid 

where we used the fact that all the baryons are contained well inside 20kpc, and 

Cl = 2-f^f~2, Al = ^~2, (45) 
a log it a log g 

Note that we do not have freedom in chosing the law of the gravity, or £i or Ai: they are all determined by the nearly 
flat circular velocity curves in the Milky Way and Andromeda, which are inferred from the fairly constant velocity 
dispersions of the outer satellites and other tracers of the outer halo such as PNe. 

Satellites with known proper motions often suggest only mildly eccentric orbits: pericentric distances are rarely 
smaller than 1/2 of the present distances, and almost never smaller than 1/5 of the present distance (Dinescu et 
al. 1999, 2004 and references therein). Outer globulars on long period orbits also have time to fill the instantaneous 
Roche lobe by two-body relaxation (Bellazzini 2004). Globular clusters as a class generally exhibit extra-tidal stars 
under deep observations; a sharp break of the star count profile and/or sometimes two-dimensional confirmation of 
unvirialised structures outside the King radius have been shown for globular clusters in the Milky Way, Andromeda, 
and other galaxies (Leon, Meylan & Combes 2000, Lehmann & Scholz 1997, Sohn et al. 2003, Harris et al. 2002). To 
be generous we assume the instantaneous Roche lobe has a filling factor F\ ~ 0.01 — 1, or F -1 / 3 ~ F 1 1//3 ~ 1 — 5. So 
equation l|44[) predicts that the majority of satellites should lie in a band with dashed line boundaries as marked in 
Fig.0 

In reality we observe a number of outliers to the inequality equation l|44|l and the expectation band. Overall we 
observe a very large scatter of O m i n for satellites in similar gravitational field with similar expected values of ©ir (cf. 
Fig. [3). Our result depends on the uncertain orbital pericentres of the halo satellites. A better understanding of the 
issue will likely come with proper motions for a large sample of halo satellites from the GAIA astrometric mission. 



4 We drop the minor factor F 2 1/3 ~ 0.8 - 1.25. 



HongSheng Zhao and LanLan Tian: Roche Lobe Shapes for testing MOND-like Modified Gravities 



11 



1.5 



c 

e 



3 



■a 
C 



0.5 




0.36 



0.38 



0.4 

Expectation 



0.42 



0.44 



Fig. 7. shows expectations vs. observations of the rescaled limiting minor axis radii m i n of outer satellites (beyond 
20kpc) of the Milky Way and Andromeda. The horizontal axis is the value of O m i n predicted from mass-tracing- 
light models assuming a satellite on a circular orbit filling its Roche lobe. The expectation band is shown by two 
slanted dashed lines; the expectation-observation equality line (thick diagonal dashed line) should be an upper limit 
for satellites are on eccentric orbits and not filling their Roche lobes. There are several outliers. To be conservative, 
we adopt a relatively low luminosity for the Milky Way and Andromeda so that = 4x ^*J L . To be accurate we 

V o 1 + ( TooT^ ) ( cf - Wilkinson & Evans 1999) for the 



adopt a slightly falling circular velocity curve V c i 
Milky Way and Andromeda. 



6. Conclusion 

The shapes of Roche lobes of a binary in a class of gravity theory including MOND are given analytically. The volume 
ratio of the Roche lobes is proven to scale linearly with the true mass ratio of the binary, hence models with dark 
matter would predict a Roche lobe volume V times that predicted in mass-tracing- light models, where the factor 
^ = (x)sat (x)h var i es from V ~ 1/10 for typical outer globular clusters in a host galaxy to V ~ 10 for typical 
dwarf galaxies in a host galaxy. Observed limiting radii of outer satellites of the Milky Way have a large scatter, 
with some outliers by a factor of 4 above or factor 10 below predictions of our basic models for the Roche lobes in 
MOND-like mass-tracing-light gravity theories. The flattenings of the Roche lobe are independent of the mass ratio, 
but are sensitive to the function n(g) in modifications to the law of gravity. Precise measurements of the exact shapes 
of other Roche-lobe-filling systems (e.g., mass-lossing stars or gas clouds in the outer halo of the Galaxy) might also 
test the law of gravity in the weak regime. 
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